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Abstract 

' 

^— ^ ' Results for light quark masses obtained from lattice QCD simulations 

, are compared and contrasted with other determinations. Relevance of 

f"^ ■ these results to estimates of e'/e is discussed. 



1 Determination of Light Quark Masses 

Quark masses are not physical observables in QCD, rather they enter the theory 
as parameters in the Lagrangian. Their values depend on the QCD renormal- 
isation scale, and three quantitative approaches have been used to determine 
them. 



1.1 Chiral perturbation theory 



This is a low energy (E < 1 GeV) effective field theory of QCD in presence of 
spontaneous chiral symmetry breaking. With m q <C A-qcd , the pseudo-Nambu- 
Goldstone bosons are the dominant fields at low energy. 0(p 2 ,m q ) terms in 
the effective Lagrangian are fixed by the spectrum. 0(p 4 ,p 2 m q ,m q [ ) terms are 
estimated using resonance saturation and large N c power counting, as well as 
from phenomenological fits to various form factors. Small electromagnetic and 
isospin breaking effects are systematically included, and bounds are obtained 
on quark mass ratios. 

The range of validity of this effective field theory expansion cannot be con- 
vincingly established. Also the absolute mass scale has to be fixed from experi- 
mental data. Yet renormalisation group invariant dimcnsionless mass ratios can 
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be tightly constrained (!]]. 

— = 0.553 ±0.043, — = 18.9 ±0.8, — = 34.4 ±3.7, (1) 

m d m d m u 

2m, m, m, — to„w , „ 

— = — - = 24.4 ±1.5, — — = 40.8 ±3. 2. (2) 

m u + md m u d ma - m u 

These mass ratios are typically combined with m s (l GeV) extracted from QCD 
sum rules to give individual quark masses. 

1.2 QCD sum rules 

In this approach, two-point hadronic current correlators are evaluated using the 
Operator Product Expansion (OPE) and perturbative QCD in the Euclidean 
region q 2 < 0. The expansion in inverse powers of — q 2 is then analytically 
continued to the physical region q 2 > 0, and matched to the experimental data 
using dispersion relations. The experimental data is organised in terms of con- 
tributions from well-known leading poles and subleading branch-cuts from the 
QCD continuum. Subtractions are used in the dispersion relations to suppress 
large continuum contributions, while threshold factors and moments help in 
suppressing contribution of regions near the poles where perturbative QCD is 
inapplicable. The results critically depend on proper choice of the boundary 
conditions. 

Pscudoscalar, scalar, e + e~— annihilation and r— decay sum rules have been 
used by various groups to extract m s . In the MS scheme, the required (3 and 7 
functions are known to 4— loop precision, giving the most accurate results 0: 

m s (l GeV) = 162.5 ± 15.5 MeV , m s (2 GeV) = 117.8 ± 12.3 MeV . (3) 

Analysis of Cabbibo-suppressed t— decay gives m s {M T ) — 119 ± 24 MeV ||. 

The major uncertainty in the sum rule approach is the asymptotic nature 
of the perturbative QCD expansion — in the longitudinal channel it cannot be 
pushed beyond O{o 3 s ). 

1.3 Lattice QCD simulations 

These provide a first principle non-perturbative determination of quark masses 
without additional assumptions. It is convenient to convert the running quark 
masses to renormalisation group invariant quark masses, completely within the 
lattice regularisation framework. 

mf G/ = lim {(2/3o ff 2 (M)r W2/3 WM)} ■ (4) 

With improvements in simulation algorithms and computers, the systematic 
errors from finite lattice spacing, finite lattice size and the quenched approxi- 
mation are gradually coming under control. Most simulations are in the range 
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m s /3 < m q < 3m s , and m u — md is assumed. Results from several simulations 
are combined to extrapolate to fi — * oo. 

Pcrturbative relations connecting lattice and MS values show that including 
higher order corrections decreases m q . Moreover, at (J, = 2 GeV, the N 3 LO 
correction is comparable to the N 2 LO one, and decreases m q by about 10% Q. 
A rigorous comparison of lattice and MS values should therefore be made at 
/i > 2 GeV, contrary to common practice. 

Here is a summary of results presented by various groups at LAT99. 

• CP-PACS: RG-improved gauge field action and tadpole-improved clover fermion 
action are used. Ward identities are used to extract m q . Nf = and Nf = 2 re- 
sults are compared, demonstrating that going to Nf = 2 reduces the systematic 
error in the lattice scale. The results converted to the MS scheme are 0: 

m ud {2 GeV) = 3.3(4) MeV , m s (2 GeV) = 84(7) MeV . (5) 

• ALPHA/UKQCD: Schrodinger functional method is used with 0(a) improved 
action and Nf = 0. PCAC relation is used to extract m q . With (M^ro) 2 = 
1.5736 fixing the reference scale, the results are ||: 

2M ref = mf GI + m£P = 143(5) MeV , m s {2 GeV) = 94(4) MeV . (6) 

• QCDSF: Schrodinger functional method is used with 0(a) improved action 
and Nf = 0. PCAC relation is used to extract m q . Using ro = 0.5 fm as the 
reference scale yields 0: 

m ud {2 GeV) = 4.4(2) MeV , m s (2 GeV) = 105(4) MeV . (7) 

• Rome group: 0(a) improved action is used with Nf — 0. Ward identity for 
the renormalised quark propagator is used to extract m q . With pseudoscalar 
and vector meson masses fixing the lattice scale, the results are ||: 

m ud (2 GeV) = 4.8(5) MeV , m s (2 GeV) = 111(9) MeV . (8) 

Numerical results show that, for a fixed value of the lattice scale, increasing 
Nf increases the coupling g 2 through vacuum polarisation. This translates into 
m q (ji) decreasing with increasing Nf. My educated guess, including the effect 
of a relatively light s— quark, is 

m s (2 GeV, N f = 2.5) = 90 ± 10 MeV . (9) 



2 Determination of e'/e 

e and e' parametrise indirect and direct CP— violation effects in neutral kaon 
decays into two pions. e arises from the mixing between CP— eigenstates Kl 
and Ks, and is experimentally measured to be 

e = 2.280(13) • 10~ 3 e 1 ^ , e w tt/4 . (10) 
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e' parametrises CP— violation in the decay amplitudes without Kl — Ks mixing, 
and in the standard model it arises from the complex phase in the CKM quark 
mixing matrix, mr— scattering data establish that the phases of e and e are 
almost identical. Recent results from NA31, E731, KTeV and NA48 experiments 
have determined (see Ref.|| for details) 

Re{e'/e) = (21.2 ± 4.6) • 1(T 4 . (11) 

2.1 Parametrisation using OPE 

All phenomenological explanations of e'/e take the CP— conserving data from 
experiments and estimate the CP— violating part. In the standard model, e 
is found by calculating the box diagram, while e' is found by calculating the 
g, 7, Z° penguin diagrams. Operator mixing and RG-evolution are crucial in the 
analysis, and electroweak contributions have become important due to large nit- 
The dominant components are the K — > mr matrix elements of the 4— fermion 
operators Qq and Qg, in the AI = 1/2 and AI = 3/2 channels. Let P^.az 
denote the ratios of the actual matrix elements to their values in the vacuum 
saturation approximation (VSA). Using commonly accepted scale parameters, 
and including isospin breaking effects with f^+jj' = 0.25(8), an approximate 
formula is || 

e - e » 13 . Im (v td v; s ) ■ (^^y ■ t^j^y 

■ [B 6 ,(i/2)(l-SW)-0.4B 8i( 3 /2) ] . (12) 

1 want to emphasise that the appearance of m s in the above formula is spurious; 
it arises because the VSA matrix element values are accompanied by quark mass 
factors in the chiral limit. It is silly to express one unknown matrix element in 
terms of two other unknowns, m s and Pi. a/, but that has become commonplace. 
To avoid unwanted systematic errors, one should therefore either (a) calculate 
the matrix elements directly, or (b) use the same calculational framework to 
evaluate both m s and the P^az- Much of the confusion in the literature has 
arisen from not following this simple guideline consistently. 

e'/e has been evaluated in three different frameworks: lattice QCD, large— N c 
approximation, and chiral quark model. Lattice QCD simulations show that 
^8,(3/2) is suppressed below 1, but provide no clean result for P6,(i/2)- With 
^6,(1/2) = 1-0(3) and Ps,(3/2) = 0-8(2), the estimates for e'/e are about a factor 

2 below its experimental value. 

2.2 Final state interactions 

Final state interactions (FSI) strongly influence K — > tttt decays, as shown by 
the experimental phase shifts. They are absent in VSA, lattice QCD, large— N c 
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approximation, and appear at subleading order in chiral perturbation theory. 
Their contribution to the decay amplitudes can be estimated using dispersion 
relation analysis of the experimentally measured phase shifts in the elastic mr 
channel. There is no doubt that incorporating the FSI boosts the AI = 1/2 
amplitude and suppresses the AI = 3/2 one, and hence increases e'/e. The 
exact value of the enhancement is debatable, since it depends on the boundary 
conditions used to analyse the experimental data. In a specific analysis ||To|| , 
the enhancement is the required factor of 2. Higher order chiral quark model 
calculations, which automatically include the FSI, also enhance e'/e close to its 
experimental value |fllf . 

In conclusion, the standard model, with a proper framework to handle non- 
perturbative QCD effects, is fully capable of explaining the observed value of 
e'/e; contributions to e'/e from new effects beyond the standard model must be 
kept < 5 • 1CT 4 . 

References 



[1] H. Leutwyler, Phys. Lett. B378 (1996) 313, [hep-ph/9602366| ; 
|hep-ph/9609467|. 



[2] S. Narison, in proceedings of "QCD99", hep-ph/9911454 



[3] A. Pich and J. Prades, in proceedings of "QCD99", [hep-ph/9909559 
[4] K.G. Chetyrkin and A. Retey, [hep-lat/9910332| . 



[5] CP-PACS collaboration, in the proceedings of "Lattice 99", 
|hep-lat/990905C| . 



[6] ALPHA/UKQCD collaborations, in the proceedings of "Lattice 99", 
|hep-lat/9909098j . 



[7] QCDSF collabor ation, 
hep-lat/990802£ . 



in the proceedings of "Lattice 99", 



D. Besirevic, V. Lubicz, G. Mar tinelli and M. Testa, in the proceedings of 
"Lattice 99" , |hep-lat/9909039 . 



[9] A. Buras, in the proceedings of "Kaon 99" , [hep-ph/9908395 
[10] E. Pallante and A. Pich, |hep-ph/9911233| . 



[11] S. Bertolini, J.O. Eeg and M. Fabbrichesi, [hep-ph/0002234 



5 



> ' 

5 



♦ RC, VWI with sea quark K = 

• RC, VWI with PQ K 
■ RC, AWI 

O qPW, VWI 
□ qPW, AWI 
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Figure 1: Continuum extrapolations of m uc t (top), and m s with m^-input 
(middle) and m^-input (bottom) [5]. Open and filled symbols are results for 
N f = and 2. 
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Figure 2: Continuum extrapolations of M re j in units of (a) ro, and (b) (Fk)r 
[6]. The largest value of a is not used in the (dashed) extrapolation to a = 0. 
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